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The Introduction. In general, the motion of charged particle in the 
external electromagnetic field is described by the full 4D Dirac equation (for 
the spin equal to 1/2). For time-independent fields Pauli obtained a famous 
beautiful nonrelativistic approximation (see P, paragraph 33). We are 
going to consider a special 2D case where the electric field is equal to zero. 
The remarkable theory of this operator starts in 1979-80. In the present work 
we found unification of the Theory of Ground States for the Pauli Operator 
with the Algebro-Geometric Spectral Theory of periodic scalar Schrodinger 
operators and corresponding theory of soliton nonlinear systems. Especially 
analog of the 2D Burgers hierarchy plays a key role here. The Aharonov- 
Bohm type terms in magnetic field like the quantized delta function also play 
important role in our construction. This work extends results published in 
[21 |3l H]. We pay special attention to the self-adjoint boundary problems 
allowing to reveal physical meaning of the Bloch functions with nonunitary 
multipliers. Such ideas were known for D = 1 since 1974, but for D > 1 
there extension never was discussed. 



History: The Ground level of magnetic Pauli operators. In 1979 
1980 three groups of authors studied the ground level using "the Factorization 
Property" of the 2D purely Magnetic Pauli Operator written in the Lorenz 
gauge + A2y = with Ai = = -i^x- (Avron-Seiler [AS], [5], 
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The Pauli Operator is a second order 2-component Schrodinger-type op- 
erator. We use a special system of units such that 

= L+^L- and -L± = {d^ + i^^y)^ + {dy - t^^^f ± A$, acting on 
the space of vector-function \1/ = 

We assume that the charge is also equal to 1. Following formal observation 
is extremely useful: 

Lemma 1 Let Q = dz — A^, = —{0^ + A^). The Scalar Operators 
are Strongly Factorized = QQ^ , L~ = Q^Q, dz = d = dx — idy,Az = 
Ai — iA2, magnetic field B = A$. 

The proof can be easily obtained by the direct substitution. 

The most interesting classes of magnetic fields are [AC] and [DN]. Fol- 
lowing results were obtained (see proof below): 

1. AC: Rapidly decreasing fields, 

\[B]\ = I J Bdxdy\ < oo. The ground states form a finite-dimensional 
space of dimension m ^ Z, 
m < [B]/27r <m + l. 

2. DN: Arbitrary periodic fields with integer- valued flux through the 
elementary cell ^ / /^^^^ Bdxdy = m E Z. 

The ground states generate an infinite dimensional subspace in the Hilbert 
Space L2(M^) isomorphic to the Landau level. The so-called " Magnetic-Bloch 
states" form a manifold isomorphic to the total space of the holomorphic 
vector bundle over the 2-torus; The projection exactly coincides with the 
so-called "magnetic quasimomentum" . The first Chern class of this bundle 
is equal to the basic cohomology class in the group (T^ , Z) , the fibre is 
isomorphic to C™. 

Remark 1 The experience of Differential Topology after Whitney, Pontrya- 
gin and Thorn led S. Novikov to the idea to study "Generic" operators in 
order to define their topological invariants. It started after this work and 
allowed to reveal topological invariants of operators in the space of quasimo- 
menta. In particular Chern numbers of the transversal dispersion relations 
were invented and studied in 1980-81 by S. Novikov and A. Lyskova-see 
13 [70| [77], [7^ . This set of ideas was partly rediscovered by physicists of the 
Thouless group few years later to explain the famous "Integral Quantum Hall 
Fenomenon" experimentally discovered at that time. 



2 



Theorem 1 In the cases AC and DN all ground states are the Instantons 
belonging to one spin-sector only: It means precisely following 

a. They satisfy to the 1st order equations Q~^ijj = in the sector for the 
case [B] > and Qip = in the sector L~ for the case [B] < 0. It is a simple 
prototype of the self-duality equation. 

b. They belong to the Hilbert Space L2(M^). 

Proof, a. The case AC. Let is a zero mode for the operator L^"^ = 
and \I/ G L2. The standard "Instanton argument" is =< QQ'^'^,'^ >=< 
Q~^'^,Q'^'^ > for the ground states. So Q"*"^ = 0. Every square-integrable 
solution of this equation \1/ G L2 defines a zero mode of the operator and 
vice versa. So our resuh follows if we will be able to find solutions from the 
Hilbert Space L2(U.'^). Let < [fi] < 00, m < ^[5] < m + 1 in the case 
AC (above). We look for \1/ as a product of two factors: one is an arbitrary 
holomorphic polynomial Pi, I = 0, m, other one depends on magnetic field: 



The growth of R obviously depends on the value of magnetic flux [B]: 



for / = 0, m — 1. So we found all ground states in the case AC. There are 
no square integrable solutions in the second sector L~ . 
b. The case DN. Solutions here also can be found in one spin-sector only. As 
above, we assume that the magnetic flux through elementary cell is a positive 
integer m G Z+. We consider rectangular lattice with periods 2uj G M, 2iu)' G 
iM. Let us define R by the formula 



where AR = B, K is elementary cell. We look for the solutions written in 
the form 



= Pi{z)exp{-R{x,y)} 
with special solution to the equation AR = B if [B] > 0: 






Q+^ = 0,^ = exp{-i?(x,y)}[e"^ JJ a{z - aj)]\ 
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Lemma 2 These solutions are the Magnetic Bloch functions, i.e. the eigen- 
vectors of magnetic translations T*,j = 1,2: 

Vector potential is chosen in the Lorentz gauge form A = {i^y, —i^x) for 
B = A$ and $ = —R as above. By definition, 

A{x + 2u, y) = A{x, y) + i{f\x, fiy), A{x, y + 2uj') = A{x, y) + i{f2x, f2y)- 

Proof. We use the standard transformation properties 

aiw + 2lu) = -e~2^("'+"V(w;), aiw + 2ilu') = -e-2*^'("'+*"'V(w;) 

where Ci'^) = <^w/o',ri = ({u),iri' = Ci^^') fo^' the rectangular lattice. Our 
theorem immediately follows from this properties by the direct calculation. 
After calculations, we are coming to the following formula for such values of 
parameters that both multipliers are unitary = \x2\ = 1: 

Rea = Re{?7i/u;[2 N^aj — I/tt / / zB{x,y)dxdy]}, 
. J Jk 

Ima = Imjr/'/cu'p ttj — l/vr / / zB{x,y)dxdy]}. 



By definition, the components of quasimomenta are pi,P2'- 

They have a form following from the formulas for the magnetic Bloch func- 
tions above: 

Pi + mix I2uj = Im(a — 77/0; aj), 

j 

P2 + mn/2u' = Re(a — rj' /u' aj). 

j 

The quasimomentum map is well-defined on the space of parameters de- 
scribing the states with unitary multipliers. It looks like 

Pi + W2 = p{ai, cbm) = -(2vrz/|A'|) ^ + const 



where \K\ = Auu' is an area of the elementary cell K. Its image is a torus 
T^, the fibre is a complex space C"*. 

The elementary arguments presented in the work [3, [8] show that this 
family describes all eigenstates of the operator nearby of the ground level 
e = 0. It follows from the properties of the elliptic operator with fixed 
unitary multipliers. Its Kernel is exactly m-dimensional because the index 
of operator with : '^'^ — \E'~ is equal to m, and the adjoint operator 
Q : — >■ \E'"'' has no Kernel. No singularities appear in this family for all 
points in the torus of quasimomenta (i.e. for all unitary multipliers). So we 
conclude that there is a finite nonzero gap A separating zero level from all 
other levels for the operator L^. 

Supersymmetry and Laplace Transformations. 

The operator = Q+ : ^+ \&" and S* : is called Super- 

Symmetry for L^. Here = 0, SL^ = L^S, S*L^ = L^S*. The "adjoint" 
supersymmetry operator is 5* = Q : — ^P"^, 
^*(^+) = 0. We have SS* + S*S = L^. 

It implies that all higher levels are 2-degenerate (the ground level is oo- 
degenerate). The Super-Symmetry operator was invented here in 1980s. 

It is in fact a partial case of the so-called Laplace Transformations known 
since XVIII Century for the scalar 2D Schrodinger operators. Every 2nd 
order operator of the (Hyperbolic) form L = d^dy + Ad^ + Bdy + V can be 
presented in the form L = {d^ + B){dy + A) + W = QQ' + W. We define 

L^L = WQ'W~^Q + IV, ^ ^ ^ = Q'^. 

We conclude that = implies = and 

L={dy + A- Wy/W){d, +B) + W={d,+ B){dy + A " Wy/W) + W 

where W = W - A^ + By + {\n W)^y. 

We are going to deal with Elliptic Case important for Quantum Mechan- 
ics, and replace dxhy d = dz and dy by d. So we have L = QQ^ + W-^ 
L = WQ+W-^Q + W and^ ^ = § with Q = d + A,Q+ = -d + A. 
For the special case of factorizable operators W = const this trans- 
formation acts on the whole spectrum of operator L. Physicists 
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identify it with "supersymmetry" for the Pauh Operator: = e\l/ 
implies Z^' = e§ for ^' = (5+\l>. 

The gauge invariants of operator L = QQ^ + W are Magnetic field 2B = 
Az — Az and Potential W. It is convenient to write Laplace transformations 
in terms of invariants only: 

W = W + i3, 13 = B + l/2A\ogW. 

It was observed that the whole infinite chain of Laplace Transformations is 
equivalent to the "2D Toda Lattice System". This observation is going back 
to the XIX Century school of Darboux (see [13] )• This school also started 
to consider the Cyclic Chains of Laplace Transformations and pointed out 
that they lead to several fundamental completely integrable systems. These 
works contains only formal calculations dedicated to the hyperbolic case (see 
historic discussion in [S]). 

In the work of S. Novikov with A. Veselov (1997) cyclic chains were con- 
sidered in the elliptic case where global restrictions play fundamental role. 
For smooth elliptic double-periodic operators it was proved in [13] that cor- 
responding Schrodinger operators L are always Algebro-Geometric, i.e. the 
collection of nonsingular complex Bloch-Floquet solutions to the equation 
L\E' = form an algebraic curve (i.e. of finite genus). However, this the- 
ory contains only Topologically Trivial operators where all coefficients are 
double periodic, and magnetic ffux through the elementary cell is equal to 
zero. In order to include topologically nontrivial operators, another type of 
Laplace chains was considered in [H]. Let Lq Li ^ ... ^ Ln he a Laplace 
chain such that the end operators Lq, Ln are both strongly factorizable up to 
constant: 

Ln = QnQt^ Lq = QQ^ + const. 

The ground level is highly degenerate here for the smooth double periodic 
magnetic field -B„. For n = 2 nontrivial cases were found such that some 
higher level is also infinitely degenerate. Probably, no more than 2 levels may 
be highly degenerate except the case of Landau operator in the homogeneous 
magnetic field. 

Algebro-Geometric (AG) Scalar Schrodinger operators. 
Question: Is Theory of Ground States for the Purely Magnetic Pauli 
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Operator related somehow to the Algebro-Geometric Theory of the scalar 
periodic 2D Schrodinger Operators based on the Selected Energy Level and 
2D Soliton Hierarchies? 

Let us remind here the history of this subject. 

History: The Algebro-Geometric Theory of the 2D Second Order Scalar 
Schrodinger Operators and Corresponding Soliton Hierarchies based on the 
selected energy level was started in 1976 by Manakov [TS] and Dubrovin, 
Krichever, Novikov [T6] . 

The magnetic field is always Topologically Trivial in this theory (i.e. coef- 
ficients of vector-potentials are double periodic). The idea to use following 
systems instead of ordinary Lax pairs was invented in [13] 

Nontrivial example with both operators L, H of order 2 was demonstrated 
showing that this new approach is nonempty. Long time this specific system 
was not studied. A "2-point scalar Baker-Akhiezer function" was invented 
in [16] for the solution of inverse problems for L based on one energy level. 
It allows to construct the whole hierarchy of higher systems associated with 
the original one. In particular, the whole hierarchy is completely different 
from the KP-type situation: it depends on 2 infinite sets of times t'j,tj as- 
sociated with 2 "infinite" points. Our Riemann surface here is a family of 
all nonsingular complex Bloch-Floquet functions of one selected energy level 
Lijj = ("the Complex Fermi Curve"). In the Algebro-Geometric Case 
it has finite genus. Two "infinite" points are selected on this algebraic Fermi 
curve. 

Some important time- invariant Reductions (i.e. invariant subsystems or 
sub-hierarchies) were actively studied in 1980s. Several authors found them 
either for Nonlinear Systems or for Inverse Spectral (Scattering) Data (or 
for both). Solution of this problem for the Inverse Data is more difficult: It 
implies in particular the description of all reduced hierarchy. However, the 
existence of time-invariant reduction is much easier to observe on the level 
of equation. 

Our Main Goal here is Quantum Mechanics and Spectral Theory. We use 
Nonlinear Systems as a tool for this goal. 

1. The Data leading to the Self- Adjoint Periodic Operators were found by 
Cherednik in 1980 [H]. 

2. The Data leading to the traditional operators L = —A + U{x,y) (with 
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zero magnetic field) were found by Veselov and Novikov in 1984 (see [TSl \W\ 

EDI EH). 

Extension of tliese results to the Rapidly Decreasing Potentials was stud- 
ied in the works of S. Manakov, P. Grinevich, R. Novikov and S. Novikov 
in the late 1980s (see [221 |23l |2l]). Krichever in |25] proved that every 2D 
smooth periodic potential can be approximated by the AG ones. 

The Problem solved recently in our work: Calculate AG Data for the 
Reduction leading to the Factorized Operators 

Solution of the Reduction Problem. 

Consider a very first Manakov's System Lf = [H, L] + fL where both 
L, H are second order operators: L = dxdy + Gdy + S, H = A + Fdy + A, 

Gt = GxX — Gyy + {F"^ / A)x — {G'^)x " + '^Sy, 

St = Syy Sxx '^{GS)x + (-^'S')y, 

with differential constraints 

Fx = 'iGy, Ay = 2Sx, f = 2Gx — Fy. 

Konopelchenko already pointed out in 1988 (see [26]) that the reduction 
S" = is time-invariant for this system. It looks like 2D analog of the famous 
Burgers system. 

We believe in the following Informal Principle: Every natural time-invariant 
reduction of completely integrable Soliton system can be effectively described 
in terms of the scattering (spectral) data and Riemann surfaces. 
How to describe this specific reduction for the Inverse Spectral Data? Is it 
possible? Making replacement x,y ^ z, z we are coming to elliptic operators 
most interesting for us. 

Our recent result [21 SI S] describes corresponding Inverse Prob- 
lem Data for the operator L with S = 0: Take Riemann Surface 
(the Complex Fermi Curve) splitted into nonsingular pieces P = P' |J P" 
with genuses g',g". They cross each other Pj = Qj, Pj G T",Qj e P', 
j = 0,...,/. (see Fig 1 for the self-adjoint elliptic case where g' = g"). 
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Take infinities ooi e F', 002 G F" with local parameters k'~^,k""^. Con- 
struct function ip = {jp'.ip") with asymptotic ip' ~ c{x,y)e'^ ^{1 + 0(/c'~^)), 
^// _ gfc 2^]^ _|_ and divisors of poles D',D" of degree + not 

crossing infinities and intersection points. In the hyperbolic case take the 
variables x, y instead of z, z. Time dynamics can be added in the standard 
way: 

^> = e'='--+E.>i{fc')^t'.(c + 0{k'-')), r = e^"'^+J^^>i('^")^*"(l + 0{k"~^)). 

The crossing property has a form 

V;'(Q,)=V^"(P,),j = 0,l,...,/. 

Theorem 2 Such Data generate a 2-point Baker- Akhiezer function = 
on the splitted surface F = F' [J F" and scalar operator L' = A + Gdz 
with S = such that L'lp' = L'^" = 0. 

The Data for Self-Adjoint Operators. For the selection of self- adjoint 
case we need to have g' = g" and to add the degenerate Cherednik-type 
restriction. One can say that it is simply a limit of Cherednik condition 
for the degenerate Riemann Surface F = F'ljF". We take infinite points 
00' G F' with local parameter k'~^ and 00" G F" with local parameter k"~^. 
An anti-holomorphic map should be given t -.T" ^ V such that r(oo") = 00', 
T*{k') = k", and r permutes the intersection points T{Pj) = Qky The inverse 
map F' — )■ F" we also call r, so we have r : F — F, = 1. The divisor of 
poles is D' G F' and D" G F" with degrees \D'\ = g' + l,\D"\ = g" where 
g' = g" and poles do not cross the intersection points and infinities. The 
total divisor of poles at the unified Riemann Surface F is D = -D' [J D" . We 
require as usually that there exists a meromorphic one-form Vt with simple 
zeroes in the points of divisors t{D) and simple poles in the infinite points 
00 1, 002. The divisor (K) consists of the zeroes of any holomorphic one- form, 
~ is a "linear equivalence" of divisors: 

D + t{D) ~ (K) + 00' + 00". 

Now let our surface be degenerate. The one-form Q also became degenerate. 
Its degeneration consists of two meromorphic one-forms u', u" in the Riemann 
surfaces F', F" with simple poles in the crossing points Pj — T{Qj) such that 
resp.u" + resQ.u' = 0. Their divisors should satisfy to the conditions: 

(c^O = D' + t{D") - 00' + ^ Q,, {uj") = D" + r{D') - 00" + P,-. 

3 j 
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After the non-unitary gauge transformation 



in r', r", we obtain a self-adjoint operator 

L = QQ^ = {d- A){-d -A), A = ^aiogc. 

Taking = L and L~ = Q^Q, we construct a Purely Magnetic Pauli 
Operator 

= QQ+^Q+Q, The Magnetic Field is real B = l/2Alogc, periodic or 
quasiperiodic and topologically trivial. It is nonsingular if c 7^ 0, so the 
operator is self-adjoint in this case. 

How to find ground states in the Hilbert space? It is very simple 
in the periodic nonsingular case c 7^ 0. Take 'ipQ = c^^^ in the first spin-sector 
L"*" because Q^ipo = 0. 

Take (^0 = c~^/^ in the second sector L~ because Q(po = 0. In the case 
of periodic smooth real c 7^ we have exactly two periodic ground 
state functions located in both sectors. They present the bottom 
of the CONTINUOUS SPECTRUM near the ground level e = 0. 
See below full description of all complex nonsingular Bloch functions of the 
ground level. Their relationship to the magnetic Bloch functions found in [7] 
in 1980 for the topologically nontrivial magnetic field also will be discussed 
below. 
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The Case of Genus zero (Fig 1) 



Fig 1 




vanishing cycles 



We take / + 1 intersection points presented as k' = kg and k" = ps in 
F', F", and divisor D' = (ai, ai) of degree / in F'. We have 



,2 wok"- + ...+wi 
{k' -a,)...{k' -ai) 



As we can see, c = wq. So c = X]l=o where Wg = PsZ — kgZ (see 



11 



Lemma 3 For every trigonometric/exponential polynomial c of the form 
above with any set of linear forms Wg and complex coefficients Ks there exists 
AG Data generating such function c. 

Let Ws = asX + (a^, G C^. 
Transformation c — c' = ce'^+'^^+i^y leads to the gauge equivalent operator 
(with the same magnetic field). 

There exist 3 types of Real Solutions: 

1. Purely Exponential Positive Case ("The Lump-type fields") 
Ks > 0,(a„/3,) e 

2. Periodic Trigonometric Real Case. It will be considered below jointly with 
the case g = 1. 

3. Mixed exponential/trigonometric case. It can be realized only if its "dom- 
inating part" belongs to the case 1. So we will not discuss it. 

The case 1. Let "the Tropical Sum" of the forms in the set {W} is 
nonnegative = maxs^as cos (p + /3s sin 0) > 0. 

Then c~^/^ is bounded in M^. 
For the angles /|^|(0) > we have a rapid decay 
^ 0, ^ oo, Let /(0) = max{J'(0), 0}. 
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Fig 2a 



c=ey+ey-2^+e-y- 




In every class c' G ce^,W' G i?^, the set of representatives c' with 
nonnegative J = J|^y,|(0) > forms a convex polytop Tc- Its inner part 
Tc C Tc consists of all c' such that I{w'} > 0. Open part is always 
nonempty for / > 2. 7^ is nonempty for / > 1. (see Fig 2b for / = 3) 
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Fig 2b 
b) R2 = (c,^ ^) 

Wj = ajX + PjH 




-W2=(0-l) 



Here + + e~^~^ = c 

Magnetic field is decaying for i? — > oo except some selected angles, it is 
a Lump Type Field analogous to the KP "Lump Potentials". A linear sum 
under the l/2Alog() reflects the Linearization of the Burgers Hierarchy in 
the variable c. 

[B]= [ [ Bdxdy = -1/2R (f I{w}{(p)d(f) + 0{R^^) 

All points in Tc define ground states in the Hilbert Space L2(M^). 
The boundary points define the bottom of continuous spectrum. 
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The Periodic Problem. 

Let lattice in be rectangular and z = x + iy. For every real periodic 
function c we can define a whole family of "possible" meromorphic Bloch 
functions 



(t{z + p + R) 



a{z + R) 

where f{z) is an arbitrary elliptic function. 

We have Q+VL,- = for L = L+ = QQ" 
tions we get Qip^^t _,_ = for L~ = Q^Q. 

Let c 7^ 0. We need only nonsingular functions, so our manifold is m G CP^ 



For anti-holomorphic func- 



r" with function ip'. 



c (or e 



Let c has an isotropic zero. We have larger family of nonsingular (or 
"weakly singular") -Bloch functions because ^/c/o^{z + R) became only 
weakly singular such that it may enter spectrum. So the full manifold is 
= CP^ X r where F is an elliptic curve. If c is a trigonometric polyno- 
mial, 1/c is not. We have c for L"*", so we have c' = 1/c for L~ . So we need 
to calculate ip' for all real periodic smooth functions c'. 



The case of genus 1. 



Fig 3 




vanishing cycles 



15 




We take elliptic curve F' = F" = C/A with Euclidean local parameters k, p 
(the point is "infinity"), periods 2uj G 'R,2iuj' G iM , n intersection points 
Qo,Qi,-,Qn e F' and Ro,...,Rn e F". Divisors D' = (Pi, P„), = P 
have degree n + 1,1 correspondingly. We have 

Here P = Pi + ... + P„, Q = <5o + ••• + Qn? sum as in C 

All singularity of the quantity c disappear after multiplication 

c = ca{z + Q + P)a{z + P). 

Take n = 1, Qo = —Qi, Ro = Qi, Ri = Qo and solution to the equation 
'^({Qo) = ViQo- We have P = Q + P in this case, so -1/2A|(t|^ = -2n6{z). 

So the Conclusion based on the case g = 1 is following: 
The least singular real Data lead finally to magnetic field B = —1/2 Ac which 
is periodic nonsingular with magnetic fiux equal to ONE QUANTUM UNIT. 
The original AG topologically trivial magnetic field B = 1/2 Ac extracted 
from our Inverse Spectral Data is always singular for = 1; it has total 
magnetic fiux equal to zero through the elementary cell and quantized S- 
singularity at the phone of nonsingular field B, located in the point P. So 
this field corresponds to the " Aharonov-Bohm" (AB) situation. 
For g > 1 the number of quantized 5-functions is equal to k > 1. We can 
get nonsingular AG operators only from genus 0. Both pieces of the 
original Riemann surface F = F" |J F' are presented in the form of /c-sheeted 
branching covering over elliptic curve F" — )■ Fq for g > 1 as it was in the 
works of Krichever dedicated to the elliptic KP (see [27]). 

The important Conclusion is following: Comparison with [7] shows 
that the Quantized 5-fiux does not affect spectrum nearby of the zero level. 

The complex Bloch-Floquet manifolds (consisting of nonsingular or weakly 
singular Bloch functions admitted by the spectral problem) for the level e = 
and genus 5^ = 1 is 

M = 
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with function 



a{z + p + R) 
a{z + R) 

Let us point out that both functions tp', ip" are originally singular here (7 = 1; 
Our next step is to reduce operator to the self-adjoint form L. For that 
we have to multiply both by the factor 1/a/c. Let us start with ip". 

The function ip'^^^l/y^ is " nonsingular" , i.e. it is singular weakly enough for 
entering spectrum in L2(M^) for the operator L with the quantized Aharonov- 
Bohm term. After singular unitary gauge transformations L — >■ L all singu- 
larity of this function disappear. So it enters the smooth spectrum of the 
operator L with magnetic field B where singularity is removed. Now about 
The function ip'/c is nonsingular but the function "(p' j ^ is singular. Af- 
ter the singular unitary gauge transformation it maps into function ip' / ^/c 
with simple pole. It satisfies to the equation 



So we conclude that the family 'ip' / of magnetic Bloch functions 
for the operator L with nonsingular topologically nontrivial mag- 
netic field has a pole in the variables {x,y). 

The function ip'/c is nonsingular. It satisfies to the equation Lq{iIj'/c) = 
in all nonsingular points where Lq = d{d + $2). At the same time we have 
$2 ~ 1/z near singularity, all other terms are smooth. So the function 
f = {d + ^z){'^' /c), c = e*, are antiholomorphic almost everywhere. How- 
ever, in the point P where singularity is located, / has a form / ~ m{k)/z. 
Therefore ip' /c ~ m{k) + zO{\) near this point. We conclude that 
Lq{%Ij' /c) = m{k)5p{x,y) because d{l/z) = {const)6. 

Reconsider now the case g = comparing it with g = 1. 
For c 7^ and g = the Bloch manifold is equal to the union F" |J F', 
and both are CP^; Let c has an isolated zero (minimum) which is isotropic. 
Magnetic field became singular, with (5-term. The extended instanton Bloch 
function ip^^f- _,_ became only weakly singular and capable to serve spectrum 
for the operator L~^. It depends on the points of the complex 2-manifold M^, 
and all possible values of the complex multipliers 

M2 ^ (xi, X2) G C* X C* 



17 



are presented in this family. The energy dispersion relation is degenerate 
here e : 0. 

We have = CP^ x Fq where Fq is an elhptic curve, 



Its singularity is in fact hke \foJa because c has a lattice of double zeroes 
like lap. Make singular gauge transformation 



(here a depends on z under the square root). This function transforms into 
the nonsingular magnetic Bloch function of the operator with topologi- 
cally nontrivial magnetic field B. Magnetic Bloch multipliers after transfor- 
mation are equal to the original standard Bloch multipliers. 

Our Conclusion is that the periodic case g = \ gives result similar to the 
special case of = where c has an isolated isotropic zero, interchanging 
sectors ± (i.e. z and z). The higher number A; > 1 of isotropic zeroes 
for = leads to the "higher rank" family of weakly singular 
instanton type Bloch functions M^'^^ . Removing (5-singularities by the 
singular gauge transformations we get smooth magnetic Bloch manifold M^^^ 
corresponding to the periodic magnetic field like in [7] with higher flux. 

We know that the algebro-geometric case simply corresponds to the case 
of trigonometric polynomials. We take rectangular lattice in the plane x, y. 
Following relation is true 



Now we choose normalization of such that M{k) = k. Let us remind that 

is "the Super-symmetry Operator" in periodic case. 

Let us extend our results to the "infinite" trigonometric series 
for g = 0. 

We use for that the formula 



We are dealing here with the new normalization: here kj are simply the 

lattice points. 

As a corollary we have 




(const(M))eP^"-^(")V(^ + u)V^/a{z). 
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Apply this result to the infinite series c — )■ c' = 1/c. It gives us a function 
for the second component L~ of the Pauli operator. 



The Boundary Problems. 

Problem: The component V of the Bloch-Floquet manifold does not 
affect the ordinary spectrum in the Hilbert space of functions in the whole 
plane M^. Can we use it for solving physically meaningful (i.e. self-adjoint) 
boundary problems? 

Let us remind following: For one-dimensional self-adjoint periodic Schro- 
dinger operators L = —d^ + u,u{x + T) = u{x), the 2-sheeted Riemann 
surface F of the Bloch-Floquet function Lip = Xip, T*ip = tplx + T) = xip, is 
such that all branching points a^, s > 0, are real. They are the boundaries 
separating the spectral zones from the forbidden zones (gaps) at the real line. 
In order to reconstruct potential, we need to know also the position of poles 
oiip{\ ,x,Xo) for A G F (normalized by condition iP{xq) = 1). They are the 
points on the Riemann surface 7j G F forming together a " Divisor D" . There 
is exactly one pole inside of every finite gap for the smooth potentials: 

Ijixo) G [a2j-i,a2j]. 

The zeroes of ip are located exactly in the points 7j(a:). Therefore our Bloch- 
Floquet function f{x*) = ip{'yj{x),x*) satisfies to the Dirichlet boundary 
problem at the interval x* G [x, x + T] where T is a period because f{x) = 
/(x + T) = 0. 

Another interpretation follows from the fact that the same function f{x*) 
decreases to the right or to the left from the zero f\x*=x = corresponding 
to the sheet of Riemann surface where the pole 7(^0) is actually located. So 
the function f{x*) belongs to the discrete spectrum on the corresponding 
half-line x* G [x, +00] or x* G [—00, x]. 

Question: Can we find similar realization for the Bloch functions 
ip', ip" for the periodic smooth topologically trivial magnetic field 
B = 1/2 A In c where c is real and nonzero? 
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The classical John von Neumann theory of self-adjoint boundary condi- 
tions (as I.Gelfand pointed out to S.Novikov many years ago, about 1971) 
follows to the scheme quite similar to the Novikov's idea to construct the 
"Hermitian K-theory" over rings with involution: Self-adjoint extensions of 
symmetric operators correspond to Lagrangian subspaces in the Hamilto- 
nian module over the ring of functions. Let 7 C be a contour in the torus 
presenting the boundary of our domain (or its connected component). 

Consider first the "scalar" boundary conditions for the operator L = 
= QQ^ or L = L~ = Q^Q separately, not mixing different components. 
We reduce the integral to the boundary 



where V'i(^) = ^ni^l-y = (c^n + A.n)'^\-t and 1IJ2 = 4'\'y Here n means the 
vector normal to the boundary component (external). To make self-adjoint 
extension means to select a "Lagrangian Subspace" {tpi{t), ip2{t)) € A for the 
skew hermitian boundary form such that 



for every pair of elements in the subspace A. There are following types of 
boundary conditions. 

The "Ultralocal" boundary conditions such that no derivatives 
along the directions tangent to the boundary are involved. They are 
elliptic. 

1. The Dirichlet b.c il)2 = ipl-y = 0. 

2. The Neuman b.c. ipi = Vnipl-f = 0. 

3. The Leontovich b.c. a(t)ilji + (3{t)ip2 = with (a, /3) real and nonzero 
for all t G 7. Only ratio {a(t)/(3{t))^^ is invariant here. This b.c. has a 
topological charge (the degree of map 7 = S*^ — t- S*^ = MP^). 

The General "Local" b.c. They are nonelliptic. 

Let a pair of scalar differential operators on the circle U, V be given 
such that one of them is invertible (for example, equal to 1). Our b.c. is 
f/^i = Vtlj2. This b.c. is self-adjoint if {UV+)+ = UV+ = VU+. In the 
special case then one of them is equal to 1, another one is an arbitrary self- 
adjoint differential operator on the boundary. 
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An interesting special case we have if f/ = 1 , and is a first order operator 
V = iadt + v(t), a, v(t) G M. Tlie b.c. of the form ipi = Vn^J = '^i'^Pt + vit)^!) 
we call the "d-bar b.c." or d-h.c. The reason for that is following: we can 
write this b.c. (for the sign — ) in the form 

[g+V' - e^^Wt;(t)^/']^ = 

where Q+ = -{dz - A^) and Ai = i^y,A2 = -i^^.A^ = is an 

angle rotating the coordinate frame x, y to the Frenet frame n, dt- The d-bar 
b.c. is self-adjoint but not ultralocal. It is nonelliptic. We could not find any 
traces of it in the literature. However, it is very closed to the elliptic b.c. and 
self-adjoint. One can imagine that such b.c. might appear in the problem 
of physics. Replacing the operator V hj V = +idt + v{t) we obtain similar 
"d— b.c." where is replaced hj Q = d — It has similar properties. 

There are also the General Nonlocal Self- Adjoint b.c. We do not 

consider here nonlocal boundary conditions. 

Boundary Conditions for Pauli operators mixing com- 
ponents. 

The operator acts on vector-functions \E' = After calcula- 

tions we obtain 

< L^vl/, $ > - < vl>, >= - ^+0^] + [i'T^2 - i^2^l]}dt 

where 

The Ultralocal b.c. They are elliptic. 

Our skew Hermitian form has signature 2, 2 in every point of boundary 
t G 7. The family of all Lagrangian subspaces is isomorphic to the group 
U2 for every t E '~f. So all collection of ultralocal self-adjoint b.c. can be 
classified by the smooth maps — )■ U2 with topological charge from the 
group vri([/2) = ^- Lagrangian subspaces in every point t G 7 may have 3 
types: 

1. Subspaces of the form 
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where matrix R is Hermitian and a, /3 = ±. 

2. Similar subspaces with numbers 1 and 2 interchanged ip2 = Ripi- 

The intersection of the types 1 and 2 is exactly a subclass such that Hermitian 
matrix R is nondegenerate. 

The types 1 and 2 can be unified by the form Uipi = Vip2 where one of 
matrices [/, V is invertible and UV^ = VU^ is Hermitian. 

3. Subspaces of the form aV'^ = 0,cip2+d4'2 — 0- Lagrangian property 
implies ac + bd = 0. We can normalize these equations reducing them to the 
form |ap + |6p = 1, |cp + = 1. We are coming to the arbitrary matrix 
from the group SU2 numerating such special Lagrangian subspaces (or to the 
arbitrary unit vector (a, 6) G determining the whole matrix). 

The General Local b.c. mixing components. 

We can write a natural differential analog of the scalar types above in the 
form 

where one of two operators U, V on the circle is invertible. Our requirement 
is UV+ = VU+. 

The most interesting for us are the following b.c where no second derivatives 
of \i/ in the boundary points are involved. 

—ip~ = cip^ + dVip' 

where A'^ = a^idr + a are the first order operators along the boundary. Here 
the functions a'^if), a(t), d{t) G -R, h{t), c{t) are smooth and c = h. 

This b.c are nonelliptic if order of operator A'^ is more than zero. The 
ultralocal case corresponds to the zero order a"*" = 0. We are going to discuss 
the following 

Problem. Let \E' = {ip^^ip') be a formal zero mode solution to the 
Pauli operator = 0. For which contours it satisfies to some self-adjoint 
boundary condition? 

We are going to consider below this Problem both in mixing and nonmix- 
ing cases. 

Our intention is to show that for some special mixing and non- 
mixing b.c. these contours can be described as a trajectories of 
dynamical system (foliation) in the 2-torus or in the covering 
space R^. Consider first the Ultralocal case. 
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Theorem 3 Let d = and = 0. The coefficients a,c^b are such that 
a ^ R and c = b if and only if the boundary curve 7 = dD is a leaf of 
foliation u = where 

CO = (iV'+l^ + \^~\'^)d^ + IV'^p * de^ + * de~ 

Here ip"^ = \ip^\exp{i6'^} and *d^ is a vector potential of magnetic field, 
star * is a standard duality operator on the complex plane. 

The case 6 = leads to the nonmixing b.c. We wiU discuss this case 
below in details. For the d-bar problem we have = —1. Assuming as 
above that d = we are coming to the similar dynamical system. It is 
investigated below in details for the nonmixing case b = c = 0. 

The proof of this theorem is not complicated. It can be obtained by the 
elementary manipulation. 

The case which is most interesting for us is t/;"*" = ip' and = Xip" 
where A is a parameter, and ip" = Q'^ip'. Here is the operator which 
is treated in quantum theory as a supersymmetry for (see above). Our 
dynamical system depends on A. For A = we have a Nonmixing Case. 

This pair defines a General Bloch Solution of the zero energy level to 
the Physical Pauli Operator with some (nonunitary) multiplier like in the 
one-dimensional case where we had Bloch solutions in the forbidden bands 
as eigenfunctions of the self-adjoint boundary problem.. 

Let us make some very general remarks concerning the Nonmixed Bound- 
ary Problems: 

Let a contour 7 C be given with following Domain D = \ 7. Find 
solution to the boundary problem 1 or 2 such that L~^iIj = in D and 

The Ultralocal Problems. Let a(t)V„'?/' + /3{t)tp = at the contour 
7. Consider first the Problem 1 with contours homotopic to zero in 
(see Fig 4 a), dD = 7. Solution of the "Leontovich type" boundary problem 
ip satisfying to this relation should be found from the integral equation 

i) = ! ! [p{k)i)' + q{k)i)"]d^k. 
J Jm? 

For the Problem 2 with contours 7 nonhomotopic to zero in we 

have dD = 7 IJ 7' (see Fig 4 c). There is a natural Z-covering D ^ D which 
is a strip in with average direction (7 = [7] G i/i(T^,Z) = 1? (see Fig 4 
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b). The element G C C = acts freely g : D ^ D as a. shift by the 
complex number g : z ^ z + g, z z + g. By definition, the components of 
boundary 

dD = 7IJ7' 

are coverings over 7,7'. g can be viewed as a vector of lattice in M^. 
Another generator g' G complementary to g maps exactly : 7 — 7'. We 
are saying that 7' (7') is located to the right from 7(7) if \x{g')\ < 1. Let 
c 7^ 0. All functions tp' (k) , tp" (k) with k orthogonal to g have unitary Bloch 
multipliers |x(5f)| = 1 along the shift g. We have 4 domains Dj^_,D_,D in 
M2 and D in T^jsee Fig 4): 

Z}+|JD_=M^ D+flD_=7, (?'P-)CD_, (7'(D+)f|D_ = I). 

We require a(t)V„^/' + P{t)il) = on the contours 7,7' but allow to have 
different pairs of coefRcients (a, {a' 13') for 7 and 7' for the domain 
D such that both of them enter the boundary. The function ip should 
be constructed essentially from ip' but we can add also ip" if necessary: The 
most general possibility is that ip is given by the integral along the /c-axis 
orthogonal to g: 

{I):iP= [ [i^'ik, X, y)p{k) + ij"{k, X, y)q{k)]dk 

Jk>0 

for the half-plane C 

(JJ) ■.tp= [tp'ik, X, y)p{k) + ij"{k, X, y)q{k)]dk 

Jk<0 

for the half-plane D_ C 

(J/J) : ^ = /" [^'(A;, X, + x, 

for the strip Z) C 

(JV^)^ : = ^[ij'{km,x,y)pm{i^) + ij"{km,x,y)]qmin) 

for D C T^, where corresponding km are such that all ip' {km) , ip" [km) have 
the same fixed unitary multiplier x = e*^™^ (i.e. km = k^ + 2711771/ g where 
m G Z). 
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The local "d-bar" (or "d-") type Boundary Problems. 

As above in the Ultralocal case we have two types of contours : 
Problem 1 with contours 7 homotopic to zero in where 7 = dD. 
Problem 2 with contours '~f (Z T"^ non-homotopic to zero in the torus; here 

aD = 7uy- 



For the second case (most interesting for us) we consider the same type 
of domains D+, D, D with boundaries non-homotopic to zero in and 
their coverings in as above (see Fig 4). 




/ 



D 



D 



g'(D_) 



Fig 4a 
7 is homotopic to 



R^ = D+UD_ D = g'(D+) n D_ 

Fig 4b 
7 is not homotopic to 

Fig 4a and Fig 4b are associated with Example 2 below. Here g is parallel 
to the average direction of 7 
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Y' 



2co 



4(0 




x=0 




x= 2 CO 



D = T'^\j, aD = 7uy 

Fig 4c Fig 4d 

Fig 4c and Fig 4d are associated with Example 1. The contour 7 is simply 
X = 0, —00 < ?/ < 00 in Fig 4c. 



Our b.c. relation along the boundary is (with possibly different coeffi- 
cients in the different components of boundary) 

or 

for real u, v. We can try to find solution ip using the integral equations as 
above. 



We already discussed above the following question: which general re- 
strictions should be satisfied for the class of functions ijj solving 
the Leontovich type problem in the magnetic field with vector- 
potential A = -i * = i{-^ydx + (^xdy)l 

We consider now only nonmixing case. Assuming that il) = p{x, ?/)e*^*^^'^^ 
we can see that the ratio Vni^/i^ should be real. Here n is a normal vector 
to the boundary contour 7 C dD. We are coming to the following statement 
(a partial case of the theorem above): 
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Lemma 4 The boundary contour 'y for the self-adjoint scalar operator Lip = 
Xip with any Leontovich type b.c. should be a leaf of foliation given by the 
equation *{d9 + A/i) = = or {9y + ^x)dx + {-O^ + ^y)dy = if the 
phase function 9{x,y) is known. Vice versa, if the contour •y is known, the 
restriction *d9 = —d^ along the contour should be satisfied. 



The Special Contours where solution to the Nonmixing B.C. is 
a Bloch function ip'. 

Example 1: ip' as a solution to the Leontovich Boundary Prob- 
lem for the Special Contours of the first kind. 

Let g = 0. We construct solution with 2n + l crossing points pj, kj where 

Pj = a -\- ibj, Pn+j — Pji kn+j = kj. 

We have constants (nothing to do with Bloch multipliers) Kn+j = Hj and all 
parameters real. Let Yl^=i S ~ 0- have 

n 

c= kq + 2 Kj cos(2(6jX + ay)), 

6 I H~ ^ I ~t~ ^'TL I j 

k \ ^ 1 ^ ^Ti+j 

For X = we have c = kq = const. We call such contours Special, of the 
first kind. We take now contour 7 : a; = or c = const with boundary 
problem 

ipn = a{x)ilj. 

Let us point out here that for such contours we have = dn because 
(i$ = at the contour $ = Inc = const and A = {i^y, —i^^). We need 
to satisfy condition ip'^/ip G M for x = 0. Presenting ip' in the form ip' = 
^kz^j: _|_ j^g-^ WioX ip'^ I ip' e M follows for the real /c G M from the 

requirements g\x=o = gx\x=Q = 0. After elementary calculations we see that 
g\x=o = Ai cos(ay) + A2 sin(ay) and gx\x=o = -Bi cos(ay) + B2 sm{ay) where 
the constants Ag, Bg are linearly expressed through the constants Kj above. 
Our boundary condition is equivalent to 4 homogeneous linear equations 
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Ai = A2 = Bi = B2 = 0. So for n > 5 we can find solution for every fixed 
values of constants a,bj,k G M. The real values of k are orthogonal to the 
direction of the contour a: = 0. So we proved following 

Lemma 5 For every set of exponents with odd number 2n + 1 of crossing 
points andn > 5, there exists a divisor (i.e. the set of constants Hj) such that 
corresponding Bloch-Floquet function if)' {k.,x,y) with real A; G M satisfies to 
the Leontovich b.c. tp'^/ijj' for the contour'-/ : x = 0. The Bloch multiplier 
along the direction of contour is unimodular n = e~^^'^. 

So our problem is self-adjoint. Corresponding function il)'{k,x,y) serves 
as a solution to the Leontovich boundary problem in the domain D: [0 < 
X < 2u] , and in one of two half-planes : [0 < x] or D_ : [x < 0] depending 
on sign of the corresponding value of k. It serves also for the spectrum in 
the domain D where x G [0, 2uj] in the torus (i.e. periodic in the variable 
y), see Fig 4c and 4d. 



Fig 5 below is associated with the next Example 2 dedicated to d-bar 
boundary conditions. It shows situation listed in the final conclusions as a 
"maximal" boundary contour homotopic to zero such that the area inside 
coincides with the whole area of the torus. In particular, the magnetic flux 
through this domain is equal to zero according to our assumptions on the 
class of magnetic fields under consideration. In this case a boundary cycle 
consists of separatrices (see Fig 5a and Fig 5b ): 

Fig 5 

y 2co' y (2co,2co') 
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Example 2: as a solution to the "d-bar" boundary problem 
for the special contours of the second kind. 

We write d-bar problem in the form Q^ij) = u{t)e''^^*''ijj where m G M and 
6{t) is an angle from the frame x, y to the Frenet frame n, dr along the bound- 
ary curve 7. Let il) = ^'{k, x,y) / y/c. We have Q~^^p/ip = e*^'=^^'^^tifc(a;, ?/) 
where Uk G M. So our contour 7 should be tangent to vector (director) field 
dr in the torus. This field depends on parameter k E C It is obtained from 
the direction dx by the rotation 6*^*=^"^'^^ in the point x,y or z = x + iy. This 
vector (director) field can be described by the zeroes of one-form 

Re[{Q+ip/ip)dz] = 0. 

Let 

n 

j 

and in the new normalization where poles are located in the crossing points 

Now we make a gauge transformation to the self-adjoint form L = QQ^, and 
%jj = ip' / y/c. We have Q~^ip/ip = —dijj'/ip'. But d{ijj' /k) = 2ce'^^ for the newly 
normalized ip'. We obtain result Q^ip/ip = —2cke^^/ip'. 

Therefore we described the set of contours for which d-bar boundary 
condition has ip'ik, x, y) as a solution. We call them the Special Contours 
of the second kind. They are leaves of the foliation on 2-torus given by 
the following equation: 

{cke'''/tp)dz + {cke^'/i})dz = 

on the contour 7. Here c is real and nonzero. 

The zeroes of tp' are singular points of this foliation. The manifold of 
zeroes Nq = {ip'{k, x, y) = 0} is compact and two-dimensional No G CxM."^ = 
{k, X, y}. The intersection of Nq with every plane k = const consists of finite 
number of isolated points. The surface No projects into the /c-plane. The 
image is a compact set A'o — )■ A^* C C. We know that A^* is compact and 
does not cover 0. For all A; G C \ N* our foliation on the 2-torus has no 
singular points. Our fiow is analytic. Such fiow is always homeomorphic 
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to the straight hne flow with some "rotation number" p{k) according to the 
famous classical results. So we have 



i/j'/ke'^'dz + ilj' /ke''~'dz = 0. 



Let us introduce a real function 




{k + l,){k-l,) 



+ a,- 



{k-lj){k + ljz) 



) 



) 



kz+kz 



Easy to check that dF = 2il)' /ke^'' and dF = 2il)' /ke^^ . So for all k e 
C \ A^* the contours are globally given by the equation dF = or F = 
const in M^. For all values of /c G C the nonsingular levels F = a define 
trajectories which do not approach singularity. All these trajectories are 
either compact and homotopic to zero or topologically equivalent in to 
the straight line trajectory around the torus with some rotation number. 
All components of the level F = const approaching singular points are either 
isolated singular points (the centers of foliation) or separatrices of the saddles 
(including degenerate saddles). So we proved the next statement: 

Lemma 6 Our foliation is equivalent to Hamiltonian system (foliation ) with 
hamiltonian F in every compact domain E in the universal covering E CM."^. 

As a corollary we conclude that 

1. All singular points are either centers or saddles (maybe degenerate). 

2. There are no limit cycles homotopic to zero in the torus, including limit 
cycles constructed from the pieces of separatrices from saddles to saddles. 

Is our foliation globally topologically equivalent to the Hamilto- 
nian foliation in the torus given by some multi-valued Hamiltonian 
(or closed one-form)? 

Can limit cycles non-homotopic to zero exist for our system? 

Let such cycle 7 be given. It may be nonsingular or consists of several pieces 
of separatrices. Anyway, the value of function F along 7 is constant, and the 
rotation number along 7 is rational. According to our assumption, a family 
of non-closed open nonsingular trajectories 71 approaches 7 asymptotically 
for t — )■ 00. We can always choose 71 such that -^(7) 7^ -^(71). 
There are two cases: 
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Case 1: The function ip' is unbounded along the closed contour 7. It is 
possible only for such contours that -^(7) = (see proof below). Any closed 
contour non-homotopic to zero in the torus with unbounded ip' presents a 
limit cycle for our foliation-see Example to the Case 1 and Fig 6 below: 

Example to the case 1: Let the functions c,F are the same as above. 
We take following parameters: 

n = 2,k = 1/2, 12 = i/2, ai = a2 = 0.2, k = 0.55i 

c = 1 + 0.4 cos(x) + 0.4 cos{y) 

Making numerical calculation, we are coming to the limit cycles 7 where 
6 = [7] G i/i(T^,Z) and 6 is a basic cycle along the y-axis (imaginary). We 
have here -^(7) = 0: 



Fig 6 




F=0 



Here the contours F = are closed and non-homotopic to 0. Bloch 
multplier in the y-direction is different from 1. Therefore differential dF 
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exponentially increases in the y-direction, and all contours F = c tend to the 
contour F = as y — )■ +00 

Case 2: The function tp' is bounded along the open contours 71. 
According to the Lemma below it is always true if all connectivity 
components of the level F = are compact in R^. The rotation number 
is rational for the closed contour 7 non-homotopic to zero. So the trajectory 
7 corresponds to the shift : — by the lattice vector with unitary 
multiplier. After the number of iterations M we get rotation number for g'^ 
equal to 1. But the form dF transforms as g*{dF) = KdF for all lattice shifts 
g with real Bloch multipliers g*il)' = xip' ., x G M. We know that F = const 
along 7i and that the form dF is invariant under the shift g'^^ with unit 
multiplier x = 1. So we see finally that F{'y) = -F(7i) for all 71 approaching 
7 if 7 is nonsingular. This is the contradiction for the nonsingular 7 because F 
is a smooth (even analytic) function. It is nonconstant along the transversal 
section to 7. It is true also for the singular 7 constructed from separatrices. 
To prove this statement we should use specific Bloch properties of function 
F written by the explicit formula, instead of dF as above. In particular, we 
have to use its behavior under the lattice translations instead of the form 
dF. 

So we conclude that our foliation cannot have nonsingular (and 
singular) limit cycles nonhomotopic to zero and such that -^(7) 7^ 0. 

Lemma 7 Let 7 be any open nonsingular trajectory such that topological 
closure of its projection 7 in the torus does not contain the levels F = 0. 
Thenip' is bounded along the contours 7. For every nonsingular periodic field 
c 7^ and \k\ big enough this condition is satisfied. 

Proof follows immediately from the formulas for the functions F and tp' 
for \k\ big enough. The function ip' looks like e'^^ multiplied by the bounded 
periodic function. So ~ eRc(fc2)_ rpj^^ function F looks like a bounded 
periodic function multiplied by the function e^^'^^'^^). So the function e^'^^'^^^ 
is bounded along the levels F = const 7^ 0. Therefore ip' is also bounded 
along this contour. For |/c| — j- 00 big enough the asymptotic of function F 
looks like {c/\k\^)e^^~^''^ , so this function is separated from zero. Lemma is 
proved. 

We are coming to the following 
Conclusion. There are only three possibilities for the leaves of this foliation 
for any A; G C: 
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1. There exists an open nonsingular trajectory 7 for which ilj'{k,x,y) is a 
solution to the self-adjoint d-bar boundary problem in the domain D G 
like cylinder, D like strip and or D_ like warped half-planes. One should 
choose only one domain or D_ where ip'{k, x, y) / a/c has decay at infinity 
transversal to the boundary (see Fig 4) 

2. There exists a "maximal" cycle 7 homotopic to zero in the torus con- 
sisting of separatrices of the saddles such that 7 bounds the whole area of 
the torus (see Fig 5). ip' gives solution to the self-adjoint boundary problem 

3. There exists a limit cycle 7 nonhomotopic to zero such that -^(7) = 
(see Fig 6). For such boundary contour ip' does not generate solution to the 
self-adjoint problem because it is has exponential growth along the boundary. 

This picture is a natural 2D generalization of one-dimensional case where 
functions f{x*) = ip{'jj{x),x*) sitting in the gaps described an additional 
Dirichlet spectrum at the interval x* G [x, x -|- T] and in the half-line W' = 
[—oo,x] or R~ = [x, +00] for the same operator L = —d^* + u{x*). 
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